Matrix Element Distributions as a Signature of Entanglement Generation 
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We explore connections between an operator's matrix element distribution and its entanglement 
generation. Operators with matrix element distributions similar to those of random matrices gen- 
erate states of high multi-partite entanglement. This occurs even when other statistical properties 
of the operators do not conincide with random matrices. Similarly, operators with some statistical 
properties of random matrices may not exhibit random matrix element distributions and will not 
produce states with high levels of multi-partite entanglement. Finally, we show that operators with 
similar matrix element distributions generate similar amounts of entanglement. 
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Entanglement, correlations between quantum systems 
beyond what is classically possible, is an essential phe- 
nomenon in quantum information processing and a neces- 
sary resource for quantum communication. In the space 
of pure states, the overwhelming majority are of high 
multi-partite entanglement with respect to a qubit archi- 
tecture Q . Such states are necessary for quantum proto- 
cols calling for random, highly entangled states including 
superdense coding [3, remote state preparation 9|, and 
data hiding schemes 

Random states can be produced on a quantum com- 
puter by applying random unitary operators to compu- 
tational basis states. However, the implementation of 
operators drawn randomly from the space of all uni- 
tary operators, the circular unitary ensemble (CUE), is 
very inefficient. Instead, other operators have been sug- 
gested as possibly efficient substitutes for the production 
of random, highly entangled states. These include quan- 
tum chaotic operators ja] and pseudo-random operators 
0,0, 111- However, these operators are not truly random 
and thus do not uniformly cover the space of pure states. 
Here we attempt to identify what statistical properties of 
random matrices lead to the production of highly entan- 
gled, random states. Identification of such a link between 
entanglement and randomness provides a deeper under- 
standing of entanglement. In addition, isolating these 
properties may focus the search for operators with the 
ability to efficiently produce random states. Such oper- 
ators need only fulfill the identified statistical properties 
and can fall short of CUE in regards to other statistical 
distributions. 

In this paper we suggest that the element distribution 
of a given operator is a vital statistical property when 
attempting to create highly entangled states. To demon- 
strate this we first show that the larger percentage of 
CUE covered by a given ensemble, the more the matrix 
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element distribution converges to CUE and so docs the 
entanglement generation. Of course, this will generally 
cause other properties to also approach CUE. We then 
isolate various statistical properties through the use of 
operators which have certain statistical properties simi- 
lar to CUE but not others. This disconnects the matrix 
elements from other statistical properties and shows their 
primacy in entanglement production. Specifically, an op- 
erator may not produce CUE-levels of entanglement if the 
matrix element distribution does not follow CUE, despite 
having other statistical similarities to CUE. Also, CUE- 
like entanglement production can be achieved with op- 
erators that have CUE-like matrix element distributions 
even if other statistical properties do not follow CUE. We 
note that not all of the operators we explore can be ef- 
ficiently implemented on a quantum computer. Rather, 
the goal is to understand how certain operators can pro- 
duce states with high levels of entanglement. This work 
concentrates on the matrix element distribution without 
exploring higher order correlations between the elements. 
These correlation may also play an important role in en- 
tanglement generation and will be the subject of further 
study. We have briefly mentioned the importance of the 
matrix elements to entanglement generation in Ref. |^. 

The operator classes used in this work to demonstrate 
all of the above, are (1) the interpolating ensembles [loj . 
a one-parameter family of ensembles which interpolate 
between diagonal matrices with uniform, independently 
distributed elements, and CUE, (2) pseudo-random oper- 
ators proposed as possibly efficient substitutes for ran- 
dom matrices, and (3) quantum chaotic operators which 
are generally known to have many statistical properties 
similar to random matrices 

CUE matrices can be generated by multiplying eigen- 
vectors of a Hermitian matrix belonging to the Gaussian 
unitary ensemble (CUE) by a random phase and using 
the resulting vectors as the CUE matrix columns 
Thus, the squared modulus, or amplitude, of CUE ma- 
trix elements follows a distribution equal to that of CUE 
eigenvector element amplitudes. Let c\ denote the A:th 
component of the /th CUE eigenvector. The distribution 
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of amplitudes, i] = |c^p, is 

pGUE{v) = iN-m-vf-', (1) 

where N is the Hilbert space dimension. In the hmit 
N ^ oo, after rescaling to unit mean, the distribution is 
given by 

PGUE{y)=e-y, (2) 

where y — Nij |l3| . Since ij is unchanged when mul- 
tiphed by a phase, the distribution, Pcue(x), of the 
rescaled amphtude of CUE matrix elements is equal 
to PauEiy)- 

As a practical measure of multi-partite entanglement 
for an n-qubit system, we explore the average bipartite 
entanglement between each qubit and the rest of the sys- 
tem [ilE3. 

9 " 

Q^2--Y,Tr[p% (3) 

i=i 

where pj is the reduced density matrix of qubit j. In this 
work, we study the distribution of Q after one iteration of 
an operator as compared to the distribution of Q for CUE 
matrices, Pcue{Q), and the average entanglement as a 
function of time, {Q{t)), compared to the CUE average 
entanglement 

{Q)cuE = iN-2)/{N+l). (4) 

CUE matrices can also be generated based on the Hur- 
witz parameterization |l2l| and a modification of this con- 
struction is used to generate the interpolating ensembles. 
The interpolating ensembles are a one-parameter, family 
of ensembles which interpolate between diagonal matrices 
with uniform, independently distributed elements, and 
CUE. The parameter, 6, can take on any value between 0, 
representing diagonal matrices, and 1, for CUE matrices. 
The exact method for the construction of CUE matrices 
via the Hurwitz parameterization and the modifications 
needed for the interpolating ensembles is reviewed in Ap- 
pendix A. Interpolating ensembles have proved useful in 
analyzing certain scattering matrices |l6lll7l | and a quan- 
tum electron pump (18| . One of the attractive features of 
these ensembles is that their eigenvalue and eigenvector 
properties have only a weak dependence on N. The ma- 
trix elements for these operators, however, show a strong 
dependence on N and this, in turn, effects their entan- 
gling power. 

Pseudo-random matrices H, 0, 13 are operators that 
were proposed as possible efficient replacements of inef- 
ficient random operators in quantum information proto- 
cols. To implement a pseudo-random operator apply m 
iterations of the n qubit gate: random SU (2) rotation to 
each qubit, then evolve the system via all nearest neigh- 
bor couplings f^. A random SU{2) rotation is described 



by Eqs. © and The nearest neighbor coupling 

operator used is: 

n-l 

C/„„c = exp(iJ^CT^(g)a^+i), (5) 

where al is the jth qubit z-direction Pauli spin opera- 
tor. The random rotations are different for each qubit 
and each iteration, but the coupling constant is always 
7r/4 to maximize entanglement generation. After the m 
iterations, a final set of random rotations is applied. 

Reference |9| discusses various statistical properties of 
both these operator classes in connection with entangle- 
ment production. Some of these properties are displayed 
here for completion. Figs, ^andl^lshow nearest neighbor 
eigenangle spacings (the more intricate number variance 
is provided in Ref. ||9|), eigenvalue element distribution, 
matrix element distribution and one-iteration entangle- 
ment distribution for the interpolating ensemble matrices 
and pseudo-random operators for constant N. As the op- 
erators cover more of CUE, S 1 and to — )■ oo, the vari- 
ous properties including the matrix element distribution 
and the entanglement generation approach CUE distribu- 
tions, as expected. We note however, that the matrix el- 
ements and entanglement generation appear to approach 
CUE more slowly than the other properties leading us to 
suspect that there may be a connection. This is born out 
by rewriting the average of Q in terms of the elements of 
the wavcfunction 

N/2 N N/2 

iQ) = <H E (ic™nc„i')-E(i^.i'K+fi'))' 

™=lri=-f+l 9=1 

(6) 

where Ci are the elements of the wavcfunction. When 
applying an operator to an intitial computational basis 
states, the output wavcfunction elements are equivalent 
to the operator elements. To properly assert the primacy 
of this connection between matrix elements and entangle- 
ment generation we investigate groups of operators that 
fulfill only some statistical properties of CUE, but not 
others, and in that way isolate the property that causes 
CUE-like entanglement generation. 

First, we show that the eigenvalue spectrum alone can- 
not be the sole cause of entanglement generation. This is 
done in two ways: by identifying a set of operators that 
have eigenvalues with statistical properties that match 
CUE but generate no entanglement, and, second, by 
identifying operators that do not have CUE eigenvalue 
properties but nevertheless generate CUE levels of en- 
tanglement. The first set is that of diagonal operators 
in which the elements are the eigenvalues of a CUE op- 
erator. When applying diagonal operators to computa- 
tional basis states no entanglement is generated. The 
eigenvector element distribution and matrix element dis- 
tribution for these diagonal operators clearly do not fol- 
low the CUE distributions. Nevertheless, the eigenvalue 
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FIG. 1: (Color online) Distributions of nearest neighbor 
eigenvalue spacings (upper-left) and eigenvector element am- 
plitudes (upper-right) for matrices of the interpolating ensem- 
bles with S = .1 (dashed), .5 (dotted) and .9 (chained). The 
N = 256 matrix element plot (lower-left) includes the distri- 
bution for 5 = .98 (light solid line). For this S the eigenvalue 
and eigenvector distributions are indistinguishable from ran- 
dom (solid line) for the resolution of the figure. The matrix 
element distribution appears to converge more slowly which 
may be manifest in the entanglement generated by operat- 
ing with 100 8-qubit 5 — .9 (O), and .98 (o) matrices on all 
computational basis states (lower right). 



spectra fulfill all statistical properties of CUE includ- 
ing nearest-neighbor spacings and higher order correla- 
tion functions. The second set of operators are cre- 
ated as follows: let I? be a unitary diagonal opera- 
tor with random phases drawn uniformly from to 27r. 
Operators U — UcueDUqjj^ have eigenvector distri- 
butions that follow CUE, matrix element distributions 
that follow CUE, and entanglement generation equal to 
CUE. Yet, the nearest-neighbor eigenvalue distribution 
follows a Poissonian and not the Wigner-Dyson distribu- 
tion ijl^ . What we see from these types of operators 
is that the eigenvalue distribution of an operator is not a 
primary factor in an operators' entanglement generation. 
This is not to say that operators with high entanglement 
generation never follow the Wigner-Dyson distribution, 
as we will see they usually do. Rather, the eigenvalue 
distribution is not the defining property generating en- 
tanglement. 

The next step is to divorce the matrix element dis- 
tribution and entanglement generation of an operator 
from its eigenvector distribution. This cannot be done 
completely. Ref. proves that the entanglement of 

the eigenvectors generally provides a lower bound for the 
asymptotic time value of entanglement generation. This 
proof is done for bipartite entanglement and can be ex- 
tended to our investigations because Q is merely the av- 
erage of bipartite entanglement between each individual 
qubit and the rest of the system. This bound, however, 
does not exclude the possibility of a high entanglement 
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FIG. 2: (Color online) Distribution of nearest neighbor eigen- 
value spacings (top left), eigenvector elements (top right), ma- 
trix elements (bottom left), and Q (bottom right) for A*' = 256 
pseudo-random maps of m = 2 (-I-), 4 (x), 8 (□), and 16 
(O)- The eigenvalue and eigenvector distributions appear to 
converge to that of CUE (solid lines) more quickly than the 
matrix element distribution and entanglement distribution. 
To approach Pcue{Q) with one iteration of a map requires 
m ~ 40 HQ. 



generation for operators without random eigenvectors. 
Nor does this result tell how long it takes to reach this 
asymptotic value, a necessary question when trying to 
create a highly entangled state on a quantum computer. 
Through the use of the interpolating ensembles discussed 
above we show how the matrix element distribution of an 
operator relates to these issues. 

As mentioned, many of the statistical properties of in- 
terpolating ensemble matrices are only weakly dependent 
on N. However, the matrix element distribution and 
the entanglement generation are strongly dependent on 
N. Fig. 131 shows the eigenvector element distributions 
and matrix element distributions for S — .9 interpolating 
operators at various values of TV. As N decreases the 
matrix element distributions approach the CUE distri- 
bution. The eigenvector elements on the other hand, are 
practically constant with N. 

The average entanglement generated in one iteration 
of the operator, compared to that expected from CUE 
matrices, is shown in Fig. 0] as a function of TV for the 
same 6 = .9 operators. As N gets smaller the average 
entanglement generated approaches that of CUE, in a 
way similar to the matrix element distribution. This de- 
spite the fact that the eigenvector element distribution 
remains constant as a function of N. For large S in- 
terpolating ensemble, say S = .99, the eigenvector and 
eigenvalue distributions will be practically indistinguish- 
able from CUE. However, as we see here, the matrix el- 
ement distribution, and thus the entanglement genera- 
tion, will fall short of CUE. The larger the Hilbert space 
the further from CUE. The interpolating ensemble op- 
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FIG. 3: (Color online) Matrix element distribution (left) and 
eigenvector element distribution (right) for interpolating en- 
semble matrices with 5 = .9 and = 256 (x), 128 (O), 64 
(+), 32 (□), 16 (o), and 8 (A). The matrix element distri- 
butions gets further and further from that expected of CUE 
(dashed line for the A*' — > oo limit and chained line for A'^ = 8) 
as A'' increases. This is in contrast with the eigenvector ele- 
ment distribution in which is remarkably stable as a function 
of N. 



erators thus provide some divergence between the eigen- 
vector distribution and the matrix element distribution. 
For small N the matrix element distribution and entan- 
glement generation are practically random though the 
eigenvector element distribution is not. For large N and 
large S the eigenvector element distribution is practically 
random but the matrix element distribution and entan- 
glement generation are not. This further demonstrates 
the importance of an operators' matrix elements in en- 
tanglement generation. The time evolution of these op- 
erators is explored below. 

A quantum computer programmer starting with a 
computational basis state and attempting to generate 
a random state of high multi-partite entanglement will 
want to apply an operator with a matrix element distri- 
bution as close as possible to CUE. Most likely the oper- 
ator will also exhibit other statistical properties close to 
CUE, as occurs with pseudo-random operators 0, but 
that need not be the case. 

Identifying one statistical property that is the domi- 
nant cause of high entanglement generation is also impor- 
tant for understanding entanglement as a quantum phe- 
nomenon and its relation to quantum chaos. Quantum 
chaotic operators are known to exhibit many properties 
of random matrices 0, H^l including the ability to pro- 
duce entanglement. Numerical simulations of two cou- 
pled subsystems demonstrate the greater entanglement 
generation of chaotic versus regular quantum dynamics 
IE EE 1^ EE and analytical results have been obtained 
through various methods EE EE EE ■ As with the previ- 
ous classes of operators, we are interested in a quantum 



FIG. 4: (Color online) Entanglement spectra for CUE op- 
erators (lines) and 5 — .9 interpolating ensemble operators 
(shapes), for N = 8 (solid line and x), 16 (dotted line and 
O), 32 (chained line and -I-) and 64 (dashed line and □). As A' 
increases the distributions of the interpolating ensemble oper- 
ators diverges further from the CUE distribution. The inset 
shows the difference between average entanglement generation 
for CUE operators and the 5 = .9 interpolating ensemble op- 
erators, AQ, as a function of number of qubits, n = log2{N). 
As N increases the average entanglement production gets fur- 
ther from the CUE average. 



chaotic operator's ability to produce highly entangled 
states from initial computational basis states. Apply- 
ing a chaotic operator once will not, in general, produce 
entanglement on par with random operators 0, |E ■ This 
is in line with the deviant short time behavior of chaotic 
systems with respect to other statistical properties such 
as the level or number variance. The deviant behavior is 
attributed to short periodic phase space orbits EE ■ Long 
time entanglement generation behavior is related to the 
operator's eigenvectors EE- We have already demon- 
strated that the operator eigenvalues and eigenvectors 
are not primary with respect to entanglement genera- 
tion, thus we demonstrate the observed short and long 
time entanglement generation behavior is also reflected 
in the matrix element distribution. 

Upon increasing the number of iterations, i, of a quan- 
tum chaotic operator the average entanglement of initial 
computational basis states, {Q{t)), can approach that of 
random operators. Similarly, the matrix element distri- 
bution of chaotic operators at higher powers approaches 
the CUE distribution. To demonstrate this we revisit 
the entanglement production of the '30'| quantum baker's 
map Pi and explore other quantized chaotic maps. 

Initial computational basis states evolved under the 
quantum baker's map attain (Q) values close to (Qcue) 
only at large t 0. This is understood based on the 
baker's map matrix element distrubtion which does not 
at all resemble Pcue{x), Fig. [Sp. However, for t = 100 
the distribution is much closer to the CUE distribution. 
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For an 8 qubit map {Q{t = 1)) is only .3080, compared 
to {Qcue) = .9883, while {Q{t = 100)) is .9597. It 
is important to note that the quantum baker's map for 
Hilbert space dimensions which are a power of 2 is known 
to have an almost Poissonian nearest neighbor eigenvalue 
spectrum [s^, |^ • Thus, at long times we see relatively 
high entanglement generation without the presence of 
a Wigner-Dyson distribution. Rather, more iterations 
lead to increased matrix element randomness causing the 
greater entanglement generation. 

We study two other examples of quan tized chaotic 
maps: the quantum sawtooth map [S^, |33 , 



Us, 



,-i7r/4 



N 



^ik-Km /N ^17^ {n—m) /N 



and the quantum Harper map |34| . 
UH^e 



_ ^iN^ cos{2Trq/N) ^iN^ cos{2Trp/N) 



(7) 



(8) 



All elements of the chaotic, k = 1.5, and regular, k — 
— 1.5, sawtooth maps have equal amplitude. One itera- 
tion of either map on any computational basis state yields 
a state with Q = 1. For the chaotic sawtooth, the matrix 
element randomness increases with t, such that at i = 50 
the matrix element distribution is practically Pcue{x) 
and {Q{t = 50)) = .98826. For the regular sawtooth 
(Q) oscillates wildly as seen in figure [3 This stems from 
the lack of an asymptotic randomness for the matrix el- 
ements. 

The matrix elements for the chaotic Harper, 7=1, 
deviate only slightly from Pcue{x), and {Q{t — 1)) = 
.9814. For t = 50 there is an increase in matrix ele- 
ment randomness and {Q{t — 50)) = .9882. The regu- 
lar Harper map, 7 = .1, matrix element distribution and 
(Q) also approach asymptotic limits as t increases. These 
limits fall short of the random matrix statistics but the 
average entanglement is still {Q(t 00)) ~ .95. Note 
that the asymptotic average entanglement of the regular 
Harper is about the same as that of the baker's map. 
In addition, the average entanglement after one iteration 
of the map is higher for the regular Harper than for the 
baker's map. This appears to be an exception to the 
conjecture that entanglement is a signature of quantum 
chaos. The quantum baker's map is widely considered 
chaotic since it is the quantum analog of a chaotic map. 
The Harper's map for k = —1.5 is not chaotic since it 
is the quantum analog of a regular map. Yet the entan- 
gling power of the Harper's map appears to be at least 
as good, if not better, than that of the baker's map. 

Finally, we return to the interpolating ensembles ma- 
trix element distribution and (Q) now as a function of 
time. As shown in [2ll | we expect that the entanglement 
will approach the entanglement of the eigenvectors. How- 
ever, we show here how matrix elements affect this and 
how operators with similar matrix element distributions 
lead to similar average entanglement generation. Fig. El 
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FIG. 5: (Color online) Average entanglement, (Q), over all 
8-qubit initial computational basis states as a function of time 
for quantum sawtooth maps, k = 1.5 (chained line) and k — 
— 1.5 (dotted line), and Harper maps, 7 = 1 (solid line) and 
7 = .1 (dashed line), compared to the random matrix average 
(horizontal dashed line). The chaotic maps quickly approach 
the random matrix average while the regular maps do not. 
The insets show matrix element distributions for the regular 
(light) and chaotic (dark) sawtooth maps at t = 50 (A), the 
regular (light) and chaotic (dark) Harper maps at t = 50 (B), 
and t — 1 (hght) and 100 (dark) of the baker's map (C). 



shows {Q{t)) for S = .9 (x), and .98 (Q) operators, al- 
ready explored in |^ , and the matrix element distribution 
at the same points in time. As the number of iterations 
increase the entanglement production exponentially ap- 
proaches the CUE value (dashed line). The matrix ele- 
ment distribution also converges to the CUE distribution. 
Comparing the entanglement produced by these opera- 
tors we note that (0(t = 5, 10, 15, 20, 30)) of the S = .9 
operators equal {Q{t = 1,2,3,4,6)), respectively, of the 
S = .98 operators. The matrix element distributions pro- 
ducing these entanglement values are practically equal. 

In conclusion, we have explored the connection be- 
tween an operators matrix element distribution and its 
multi-partite entangling power on initial computational 
basis states. We have shown that operators with CUE 
distributions of eigenvalues and eigenvectors are not the 
sole cause of CUE-like entangling power. CUE-like en- 
tangling power cannot be achieved without a CUE distri- 
bution of matrix elements. In addition, operators with- 
out CUE distributions of eigenvalues and eigenvectors 
can still have CUE-like entangling power if the operators 
have a random matrix element distribution. It also ap- 
pears that operators with similar matrix clement distri- 
butions generate similar amounts of entanglement. This 
analysis should provide a more specified goal in the search 
for efficient means of random state production. 

The authors thank K. Zyczkowski for clarifying in- 
terpolating ensemble generation and helpful discussions. 
The authors acknowledge support from the DARPA 
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FIG. 6: (Color online) Average entanglement generation as a 
function of time for 50 interpolating ensemble operators with 
6 = .9 (x), and .98 (O) operators on initial computational 
basis states. For both sets of operators the average entangle- 
ment produced approaches the CUE average (dashed line) as 
an exponential. The insets show the matrix element distri- 
butions for the same set of operators as a function of time. 
The left inset shows the distribution for 5 = .9 operators for 
t = 1(0), 5 (□), 10 (x), 15 {+), 20 (A), and 30 (o). The 
right inset shows the matrix element distributions for 5 — .98 
operators and f = 1(0), 2 (□), 3 (x), 4 {+), 5 (A), and 6 (o). 
Both sets of operators matrix element distributions approach 
CUE as time increases. 
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APPENDIX A 

CUE construction based on the Hurwitz parameter- 
ization starts with elementary unitary transformations, 
_£(*'■')((/),■!/>, x), with non-zero elements ^3 



kk 



= 1, k = l,...,N, k^i,j 



= e^^cosc 



smc 



E. 



S^;-^) =e-^'^cos( 

J J 



(9) 



which are used to form — 1 composite rotations 

El - i?(^-^^^'(0oi>oi,Xi) 

E, = i?(^-2'^-l)((Al2,^12,0)i?(^-l-^)(<^02,V'02,X2) 

En-1 = E^^-'^'^{(t)N-2,N-l,tpN^2,N-l,0) X 
E'-^-^\(j)N^3,N^l,1pN^3M-uO) X 

... i?(^-i^^'(0o,iv-i,^o,Jv^i,Xiv-i). (10) 

A CUE matrix is finally attained by Ucue = 
e'^"EiE2 ■ ■ ■ Em-1- The Euler angles -0, x, and a are 



drawn uniformly from the intervals 

< -0^5 < 27r, < < 27r, < a < 27r, (11) 
and (j)rs = sin^"'^(^rs"'^^'"^'^^^'')j with ^^s drawn uniformly 



from to 1. The 2x2 block s'inn with m,n 



and r = is a random SU(2) rotation with respect to 
the Haar measure. The interpolating ensembles [10| fol- 
low the same construction but the angles are drawn from 
constricted intervals 



0<^rs < "^ttS, 0<Xs < "^ttS, <a< 2nS, 

(12) 

with (j)rs = sin ^(iJ^rs^^'^'^^^'') and drawn from to 1. 
The whole is multiplied by a diagonal matrix of random 
phases drawn uniformly from to 2tt. The parameter S 
ranges from to 1 and provides a smooth transition of 
certain statistical properties between the diagonal circu- 
lar Poisson ensemble (CPE) and CUE [l^. 
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